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Abstract 

Blast  waves  (BWs)  form  when  the  wave  speed  of  an 
initially  diffusive,  supersonic  radiation  wave  becomes 
subsonic  and  creates  a  radiographically-visible, 
hydrodynamic  shock  wave.  BWs  are  a  novel  diagnostic  in 
radiation-flow,  code  validation  experiments  that  use 
Sandia’ s  Z-accelerator’s  dynamic  hohlraum  (DH)  as  a 
radiative  source.  The  physics  models  being  tested  are 
sensitive  to  delivered  energy  and  power  changes  of  better 
than  ±10%;  therefore,  precise  in-situ  radiative  power  and 
energy  measurements  are  required  for  quantitative 
comparisons  between  simulation  and  experiment.  The 
energy  sensitive  BW  diagnostic  complements  bolometric 
and  x-ray  radiometric  diagnostics  in  providing  these 
measurements.  Recent  comparisons  between  BW 
qualification  experiments  and  simulations  have  revealed  a 
spatial  dependence  on  the  radiation  source.  We  discuss 
the  experimental  design  and  sensitivities  for  the  BW 
diagnostic  and  experimental  results  in  comparison  to 
simulations  and  other  diagnostics. 

I.  INTRODUCTION 

Blast  waves  are  commonly  observed  in  nature  as 
associated  with  supemovae  explosions,  coronal  mass 
ejections,  and  comet  or  meteor  impacts  on  planets.  As  an 
example  for  supemovae,  it  is  hypothesized  that  this 
phenomenon  occurs  as  the  shock  break  outs  into  the  low 
density  gas  region  surrounding  the  imploded  star  [1-2]. 
There  the  radiation  emitted  by  the  hydrodynamic  shock 
becomes  supersonic  and  out-mns  the  shock  front  due  to 
the  lower  sound  speed  and  opacity  in  that  region  [2].  As 
the  radiation  wave  expands,  its  energy  density  lowers  to  a 
level  where  it  travels  transonically  and  a  hydrodynamic 
wave  can  be  seen  as  a  density  perturbation.  This  radiation 
driven  blast  wave  can  often  be  seen  as  the  emission 
increases  with  density.  Recent  developments  in  pulsed 
radiation  sources,  such  as  those  found  at  high  fluence 
laser  facilities  and  high  power  z-pinches,  have  enabled 
controlled  measurements  of  blast  wave  properties. 


Models  for  radiation  transport  through  diffusive  media 
are  dependent  upon  the  speed  of  radiation  relative  to  the 
hydrodynamic  sound  speed  of  the  heated  material  [3].  The 
simplest  case  is  when  the  radiation  wave  is  supersonic. 
The  radiation  wave  propagates  and  heats  new  material  too 
quickly  to  result  in  significant  material  motion.  The 
pressure  discontinuity  at  the  interface  between  the  heated 
and  unheated  material  moves  before  the  material  can 
compensate  by  expanding  into  the  unheated  region. 

Using  a  simple  one-dimensional  power  law  scaling,  the 
radiation  wave  speed  Urad  can  be  estimated  as 

_\9/ 

Urad  ^  ^  ’  where  s  is  the  linear  energy  density  of 

the  wave  and  t  is  the  time  [3].  The  material’s 
hydrodynamic  sound  speed  can  also  be  estimated  by 

cc  t  In  the  case  of  unconstrained  radiation  flow, 

the  diffusive  radiation  wave  slows  with  decreasing  energy 
density  due  to  its  spherical  expansion  and  losses  to 
absorption  into  the  medium.  When  the  sound  speed 
becomes  comparable  to  the  radiation  wave  speed,  i.e.,  Urad 
-  Us,  a  radiation-driven  blast  wave  or  Marshak  wave 
forms,  see  Figure  1  [4].  The  material  near  the  front  of  the 
radiation  wave  begins  to  move  into  the  colder  medium  in 
front.  The  hot  material  builds  up  a  density  perturbation  at 
the  pressure  discontinuity. 


Time  (AU) 

Figure  1  A  notional  plot  highlighting  how  a  radiation 
wave  forms  a  blast  wave  in  a  material. 
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Figure  2  A  silica  aerogel  foam  rests  above  one  axial 
end  of  the  dynamic  hohlraum. 

This  model  simply  illustrates  the  principal  of  the 
radiation-driven  blast  wave  phenomenon,  and  the  reader 
is  encouraged  to  read  more  extensive  derivations  [2-5]. 
The  analytic  models  are  not  easily  solved  through  the 
transonic  regime  where  the  radiation  transitions  from 
supersonic  to  subsonic  radiation  transport. 

II.  EXPERIMENT  SETUP 

These  experiments  are  performed  at  the  11.5  MJ  Sandia 
National  Laboratories-  New  Mexico  Z-machine  pulsed 
power  facility  [6].  The  experiments  use  the  Z  dynamic 
hohlraum  (DH)  as  an  energetic  x-ray  source  [7].  This 
source  starts  by  passing  current  through  two  concentric 
(nested)  arrays  of  tungsten  wires.  The  outer  array  at  20- 
mm  radius  uses  two-hundred  forty  7. 5 -pm  diameter  wires, 
while  the  inner  array  at  10-mm  radius  uses  one  hundred 
twenty  7. 5 -pm  diameter  wires.  These  tungsten  wires  are 
weighted  from  above  to  retain  tension  on  the  arrays.  The 
system  has  a  return  current  conductor  at  2 5 -mm  radius. 
These  wires  vaporize  and  implode  cylindrically  upon  a 
9. 5 -mm  high,  5 -mm  diameter,  14  mg/cm^  poly  4-methyl- 
1 -pentane  (TPX  or  ~  CH2)  plastic  foam.  When  the 
tungsten  hits  the  foam,  a  cylindrically-convergent 
radiating  shock  is  generated  that  emits  70-120  kJ  of  quasi- 
blackbody  x  rays  in  each  axial  direction  in  a  -10-ns  pulse 
with  peak  powers  of  1 1-15  TW. 

The  blast  wave  targets  rest  above  a  4-mm  diameter 
radiation  exit  hole  on  the  axial  end  of  the  TPX  foam,  see 
Figures  2  and  3.  As  radiation  flows  upward,  it  passes 
though  a  1.7-mm  high,  tapered,  25-pm  thick  gold  wall 
cone  that  is  filled  20  ±  3  mg/cw?  silica  aerogel  (Si02). 
Above  this  cone  is  a  20  ±  3  mg/cm^  silica  aerogel  filled, 
1-mm  high,  2. 4-mm  inner  diameter,  25-pm  thick  gold 
wall  cylinder.  On  the  cylinder  rests  a  4-mm  diameter  gold 
platform  that  supports  a  higher  density  (40-60  mg/cm^) 
silica  aerogel.  This  aerogel  is  the  region  where  the  blast 
wave  forms  after  radiation  emerges  from  the  2. 4-mm 
diameter  cylinder.  A  spatial  fiducial  grid  is  located  on  the 
target  to  provide  a  spatial  calibration  for  radiography. 


Blast  wave  target's  silica  aerogel 


Dynamic  hohlraum 


Figure  3  Photograph  of  the  blast  wave  target  above 
the  dynamic  hohlraum. 

A  5 -channel  x-ray  diode  (XRD)  array  views  the 
downward  propagating  axial  emission  from  the  dynamic 
hohlraum  source  [8].  This  indirect  measurement  of  input 
radiation  has  significant  variability  due  to,  possibly,  the 
influx  of  tungsten  plasma  across  the  viewing  aperture  [7], 
current  contacts  between  the  wires  and  the  transmission 
line  [7],  or  a  strong  radial  electric  field  [8].  The  measured 
temperatures  and  time-integrated  energy  depositions  are 
inconsistent  with  observations  of  target  dynamics  above 
the  dynamic  hohlraum.  In  follow-on  experiments,  a  1-mm 
diameter  hole  in  the  taper  region  of  the  target  permits  a 
view  of  the  radiation  source  by  a  5  to  8  channel  XRD 
array  above  at  30°  from  vertical  [9-10].  In  general,  the 
bottom  XRD  array  peak  radiation  temperature  Trad 
measurements  using  a  blackbody  fit  ranged  from  150  to 
180  eV  with  total  energy  delivery  of  50-90  kJ.  The 
blackbody  fits  converge  in  three  to  six  iterations  with  a 
constraint  of  5%  residual  on  the  fit.  When  measured  from 
above  through  the  taper  hole  in  later  experiments  using 
the  30°  XRD  arrays,  peak  Trad  ranged  from  180  to  220  eV. 

The  primary  diagnostic  for  these  experiments  is  the  x- 
ray  crystal  imaging  system  pioneered  by  D.  Sinars,  et  ah, 
see  Figure  4  [11].  This  imaging  diagnostic  uses  the  Z- 
Beamlet  laser  [12]  with  approximately  1  kJ  of  laser 
energy  in  a  1-ns  pulse,  preceded  by  a  -0.3  kJ  prepulse.  By 
focusing  to  a  -0.1 5 -mm  diameter  spot  on  a  manganese 
foil,  helium-like  soft  x-ray  emission  is  produced.  X  rays 
with  energies  near  6.151  keV  (k  =  2.016  A)  are  focused 


Figure  4  Photograph  of  the  blast  wave  target  with 
the  x-ray  imaging  system. 
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Figure  5  The  radiographs  are  processed  to  extract  the 
edges  of  the  blast  waves. 

and  Bragg-reflected  by  a  spherically-bent  (R  =  250-mm) 
2243  quartz  crystal.  The  crystal  provides  a  magnification 
near  six  times  for  spatial  resolutions  between  10  and  20 
microns  on  a  film  pack.  The  x  rays  that  are  transported  to 
the  film  are  filtered  by  a  layer  of  aluminized  Kapton® 
poly(4,4'-oxydiphenylene-pyromellitimide)  [C16H14O3] 
and  50- pm  beryllium  [13].  A  layer  of  double-emulsion 
Biomax-MS®  film  is  used  as  a  long-term  recording  media 
behind  the  beryllium  foil  [14].  Additionally,  a  Fuji  TR 
image  plate  is  placed  between  two  50-pm  thick  beryllium 
foils  behind  the  Biomax  film.  The  film  and  image  plate 
combination  provides  good  measures  of  spatial  resolution, 
noise,  and  absolute  flux.  This  Biomax  films  were  scanned 
using  a  10. 5 -micron  by  10. 5 -micron  aperture,  while  the 
image  plate  is  scanned  with  a  25-micron  aperture. 

The  x-ray  image  timing  is  given  by  measuring  the  time 
of  leakage  laser  light  that  is  transmitted  through  turning 
mirror  prior  to  the  final  optics  assembly,  and  referencing 
that  with  respect  to  the  timing  of  peak  DH  emission.  The 
radial  emission  is  measured  by  a  Kimfoil  filtered  channel 
on  a  radially-located  5 -channel  XRD  array  at  the  line-of- 
sight  (LOS)  5-6  position.  The  axial  emission  timing  is 
also  referenced  with  respect  to  the  LOS  5-6  XRD  Kimfoil 
channel.  Most  images  are  taken  5-10  nanoseconds  after 
peak  axial  power. 

III.  IMAGE  ANALYSIS 

The  images  obtained  through  the  radiography  system 
are  post-processed  using  ENVI+IDL  [15].  Preshot  images 
of  the  undriven  target  provide  spatial  calibrations,  “flat 
fields,”  and  target  anomaly  information.  Source  intensity 
variations  and  imaging  imperfections  introduced  by  the  x- 
ray  crystal  system  are  characterized  with  the  preshot 
images.  The  weighted  tungsten  wires  from  the  DH  are  in 


the  field  of  view.  These  create  poorly  focused  lines  on  the 
radiograph.  The  dynamic  images  are  then  median  filtered. 
The  Z-machine’s  x-ray  background,  obtained  on  a 
dynamic  shot  without  the  Z-beamlet  laser  driven  source, 
is  scaled  to  match  intensities  outside  the  typical  laser- 
backlighter  exposure  region  and  subtracted  from  the 
dynamic  image.  In  ENVI,  occurrence  measures  are  taken 
to  produce  a  set  of  statistical  images,  including  data  range, 
mean,  variance,  entropy,  and  skewness.  These  statistical 
images  are  then  recombined  into  a  multi-channel 
synthesized  multicolor  image.  A  spectral  angle  mapper  is 
then  used  to  obtain  texture  correlations  between  image 
features  and  to  generate  a  new  enhanced  image,  see 
Figure  5.  This  technique  works  extremely  well  for  low  to 
moderate  noise  images,  but  care  must  be  taken  to  retain 
the  spatial  resolution  of  the  images.  The  enhanced  edges 
are  then  extracted  for  comparison  to  simulations. 

IV.  SIMULATIONS 

Simulations  of  the  blast  wave  experiments  are 
performed  in  two-dimensional  R-Z  geometry  using  the 
radiation-hydrodynamic  code  LASNEX  [16].  The 
complicated  dynamics  of  the  DH  are  not  modeled  with  the 
simulations.  Due  to  the  DH  source’s  complexity  and 
limited  computing  capacity,  the  problem  is  run  with  a 
simple  line  source  placed  just  below  the  taper  region  of 
the  target.  The  source  is  based  upon  time  and  space 
resolved  measurements  of  the  DH  output  performed  from 
near  on-axis  using  an  x-ray  imaging  diagnostic  [17].  Both 
the  time  and  space  dependence  are  included  as  a  source  to 
the  model  [18].  Detailed  measurements  of  the  angular 
dependence  of  the  DH  source  emission  have  not  been 
made  and,  therefore,  are  neglected.  The  metrologized 
target  geometry,  material  compositions,  and  material 
densities  are  closely  replicated  in  the  R-Z  geometry  of  the 
simulation’s  initial  conditions.  Simulated  radiographs 
with  photons  at  energy  6.151  keV  are  generated  by 
postprocessing  the  simulations.  In  addition,  plots  of 
density  are  generated  to  monitor  the  position  of  the 
density  perturbation,  i.e.  the  blast  wave  position.  For 
consistency,  the  simulated  images  are  postprocessed  using 
a  similar  methods  as  those  used  to  extract  the  data  from 
experimental  images. 

V.  DISCUSSION 

We  performed  six  experiments  from  2004  through  2006 
with  40  mg/cm^  and  60  mg/cm^  silica  aerogel  foams. 
Experiments  using  40  mg/cm^  silica  aerogel  produced 
poor  contrast  on  film  for  the  blast  wave  positions.  Except 
in  one  case  where  the  wave  formed  partially  outside  the 
field  of  view,  these  images  could  be  analyzed  to  extract 
the  blast  wave  position.  Target  that  used  60  mg/cm^  silica 
aerogel  exhibited  improved  contrast  and  prevented  the 
appearance  of  blast  waves  outside  the  field  of  view.  The 
contrast  improvement  is  due  to  increased  line-integrated 
density  pR  and,  thus,  total  x-ray  absorption  in  the  shock 
front.  Also,  the  density  increase  results  in  a  shorter 
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Figure  6  Simulated  blast  wave  positions  compared  to 
experimental  measurements. 

Rosseland  mean-free-path,  which  shortens  the  distance  of 
propagation  for  the  radiation  wave.  In  all  images,  a  jet 
begins  to  appear  on-axis  far  behind  the  blast  wave.  The 
formation  of  this  jet  is  currently  being  studied. 

Spatially-dependent  source  simulations  agree  with  the 
expected  peak  source  powers  of  14-16  TW,  see  Figure  6, 
measured  using  the  XRD  array  located  at  30°  from 
vertical.  These  measurements  viewed  through  a  hole  in 
the  target’s  taper  region  on  shots  that  immediately 
followed  the  blast  wave  experiments.  In  contrast,  when  a 
spatially-uniform  source  was  used,  the  peak  power 
required  to  produce  the  blast  wave  at  the  observed 
positions  was  between  19  and  21  TW.  Such  a  power  level 
would  be  impressive,  but  is  certainly  not  realistic.  It  is 
interesting  to  note  that  the  energy  that  is  delivered  to  the 
foam  is  the  same  for  both  the  spatially-uniform  and 
spatially-dependent  cases. 

VI.  SUMMARY 

We  developed  a  blast  wave  diagnostic  using  several 
series  of  experiments  at  the  Sandia  Z-facility.  Sensitivities 
of  the  diagnostic  to  foam  densities  and  computational 
techniques  have  been  examined  in  detail.  The  blast  wave 
diagnostic  produces  measurements  in  agreement  with  the 
energy  sensitivity  predictions  reported  by  R.  Peterson,  et 
al.  in  Ref  [3]. 
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